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ABSTRACT 

The Small Magellanic Cloud (SMC) has surprisingly strong submm and mm- 
wavelength emission that is inconsistent with standard dust models, including 
those with emission from spinning dust. Here we show that the emission from 
the SMC may be understood if the interstellar dust mixture includes magnetic 
nanoparticles, emitting magnetic dipole radiation resulting from thermal fluctua- 
tions in the magnetization. The magnetic grains can be metallic iron, magnetite 
Fes04, or maghemite 7-Fe203. The required mass of iron is consistent with ele- 
mental abundance constraints. The magnetic dipole emission is predicted to be 
polarized orthogonally to the normal electric dipole radiation if the nanoparticles 
are inclusions in larger grains. We speculate that other low-metallicity galaxies 
may also have a large fraction of the interstellar Fe in magnetic materials. 

Subject headings: dust, extinction; infrared: ISM; radio continuum: ISM 



Introduction 



Low-metallicity dwarf gal axies often exhibit surprisingly strong emission at submillime 



ter and mm wavelengths (e.g.. iGalliano et al 



2010; O'Halloran et al. 2010: Galametz et al. 



2003|, 120051 : iGalametz et al.ll2009l : iGrossi et al. 



201lh . substantially exceeding what is expected 



based on the observed emission from dust at shorter wavelengths. This "submm excess" could 
in principle be due to a large mass of cold dust, but in some cases the implied dust masses 
are too large to be consistent with the observed gas mass and metallicity. 

The Small Magellanic Cloud (SMC) is a prime example of this phenomenon. The dust 
spectral energ y distribut i on (SE P) has been measured from near - infrare d through cm wave- 
lengths. Both iBot et al.l (120101 ) and iPlanck Collaboration et al.l ( 12011b! ) conclude that con- 
ventional dust models cannot account for the observed 3 mm - 600/im (100 GHz - 500 GHz) 
emission without invoking unphysically large amounts of very cold dust. 
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Large submm excesses have also been rep orted for other low-m e tallicity dwarf galaxies 



NGC 1705 has receiv ed particular attention flGalametz et al.l 120091 ; lO'Halloran et al.l 12010 



Galametz et al.ll201lf ) and substantial excesses have als o been reported for a n umber of other 
systems, including Haro 11, II Zw 40, and NGC 7674 flGalametz et al.lhoill ). 



This excess emission challenges our understanding of interstellar dust. If the submm 
excess in low-metallicity dwarfs is due to thermal emission from dust, these galaxies either 
contain surprisingly large masses of very cold dust, or the dust opacity at submm frequencies 
must substantially exceed that of the dust in normal-metallicity galaxies, such as the Milky 
Way. 

In the G alactic ISM, typically 90% or more of the Fe is missing from the gas phase 
( |Jenkinsll2009l ). locked up in solid grains. Thus Fe accounts for ~25% of the dust mass in dif- 
fuse interstellar regions, although as yet we know little about the nature of the Fe-containing 
material. Interstell ar dust models based on amorphous silicate and carbonaceous material 



(e.g., 


Mathis et a. 


Draine & Li 


2007 



1977 



Draine fe Leelll984l ; IWeingartner fe Draindl200ll ; IZubko et al.ll2004 



Draine fe Fraissd 120091 ) often posit that the Fe missing from the gas is 
incorporated in amorphous silicate material, but it is entirely possible for much or most of 
the solid-phase Fe to be in the form of metallic Fe or certain Fe oxides, such as magnetite, 
that are spontaneously magnetized. 



Draine fe Lazarianl ( 119991 . hereafter DL99) noted that fe rromagnetic or fer r imagn etic 



materials can have large opacities at microwave frequencies. iDraine fe Hensleyl (120121 ) re- 
cently re-estimated the absorption cross sections for nanoparticles of ferromagnetic or ferri- 
magnetic materials. They considered three naturally-occuring magnetic materials - metallic 
iron, magnetite, and maghemite - and found that the magnetic response implies a large 
opacity at submm and mm wavelengths. 

In this Letter we propose that magnetic nanoparticles may provide the 100-500 GHz 
opacity needed to account for the strong submm-microwave emission from the SMC. Upper 
limits on dust masses in the SMC are obtained in Section [2j The observed SED of the SMC, 
and the emission attributed to dust, is reviewed in Section [31 and in Section H] we show that 
models with Milky Way dust opacities cannot reproduce the observed SED. The contribution 
of spinning dust is discussed in Section [51 where we show that spinning dust cannot account 
for the observed emission near ~ 100 GHz. In Section [6] we consider dust models for the 
SMC that include maghemite, magnetite, and metallic iron grains. We find that the submm 
and mm excess in the SMC can be accounted for by a population of magnetic nanoparticles. 

In Section [7| we discuss other evidence for the formation of Fe or Fe-oxide nanoparticles, 
and speculate on why the dust in low-metallicity galaxies such as the SMC differs from 
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the dust in normal metallicity spirals, such as the Galaxy. We also discuss the predicted 
polarization of the emission from the SMC. Our results are summarized in Section |HJ 



Mass of the ISM in the SMC 



At a distance D = 62kpc (jSzewczyk et al.ll2009l ). the Small Magellanic Cloud (SMC) 
provides an opportunity to study the dust in a low-metallicity dwarf galaxy. The present 
study will concentrate on the 2.38° radius (f2 = 0.0054 2 sr = 17.8 deg 2 ) region (centered o n 
a 2000 = 00 h 53 m 59.6 s , 5 200 o = -72°40'16.1") studied bv lPlanck Collaboration et all fepllbh . 
The 21 cm line flux from this region at SMC radial velocities is 2100 Jy MHz (mean line in- 
tensity 1318 K km s _1 over the aperture; J. -P. Bernard 2012, priv ate communication) corre- 
sponding to optically-thin emission from M(H I) = 3.99 x 1O 8 M . IStanimirovic et al.l (119991 ) 



estimated that correction for self-absorption would raise the H I mass by 0.42 x 10 8 M(n (for 
D = 6 2 kpc). Thus we estimate M(H I) = 4.41 x 1O 8 M in the 0.00542 sr region. iLerov et al. 
koOlh find M(H 2 ) « 0.32 x 1O 8 M . Thus we take M H « 4.73 x 1O 8 M (not including He) 



within the 0.00542 sr aperture. 



Elemental abundances in th e SMC are uncertain. iRussell fe Dopital (ll992f ) estimated 
(Fe/H) SMC = 0. 2 5(Fe/ H) o , while kurt fc Dufourl f flQQsh estimated (0/H) SMC = O.2(O/H) . 



Rolleston et al.l ( 120031 ) measured the abundanc es in a main-se quence B star, and found 
(O/H) = 0.6(O/H) Q , and (Fe/H) = 0.3(Fe/H) Q . iLee et al.l fj2005h measured abundances in 3 
B-type supergiants in the SMC wing, finding (Mg/H) « O.l(Mg/H) , (Si/H) « O.2(Si/H) . 
Here we adopt an overall metallicity Z^uc ~ O.25Z . 

An upper limit on the dust mass in the SMC is obtained from the observed gas mass 
combined with the estimated abundances of elements that could form dust grains. The 
gas in diffuse HI and H 2 in the local interstellar medium is routinely strongly depleted in 
elements such as Fe and Si, with the missing material presumed to be in the form of dust. 
An inventory of the well-studied sightline toward COph allows one to estimate the dust/gas 
mass ratio bas ed on the am ount of C, O, Mg, Si, Fe, and other elements that are missing 
from the gas (IDrainel l201ll . Table 23.1)J±| If we assum e the relative abund ances of heavy 
elements in a galaxy to be similar to solar composition (jAsphmd et al.ll2009l ). then 



M dust /M H < O.OO91(Z/Z ) 
M Fe ,dust /M H < O.OO196(Z/Z ) 



(1) 
(2) 



1 Thc sightline toward COph is ice- free. We assume that ices are a negligible fraction of the total dust 
mass in the SMC. 
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where MF e ,dust is the mass of Fe contained in solid material, including ferromagnesian silicates, 
Fe oxides, and metallic Fe. The mass of Fe in magnetic materials obviously is limited by 
Mpe dust- The resulting upper limits on M^ nst and Fe in magnetic form are given in the first 
line of Table dJ 



Table 1: Dust Masses in the SMC 







Total Dust 

(10 5 M Q ) 


Magnetic Fe 
(10 5 M Q ) 


comment 


Abundance limit 




<10.7 


<2.3 




Model 1: DL07 dust, 


U min > 0.2 


13. 




violates abundance limit; poor fit 


Model 2: DL07 dust, 


U min > 0.5 


9.7 




very poor fit 


Model 3: DL07 dust 


+ 40 K Fe 


8.3 


1.4 


OK 


Model 4: DL07 dust 


+ 20 K Fe 


10.2 


2.2 


OK 


Model 5: DL07 dust 


+ 17 K 7 -Fe 2 3 


9.4 


2.2 


OK 


Model 6: DL07 dust 


+ 17 K Fe 3 4 


7.2 


2.2 


OK 



a Fe mass in magnetic material 

b for M H = 4.7 x 10 8 M Q and Z = 0.25Z ( 

a T = 40 K Fe particles 

DL07 dust with cold component 
c T = 40 K Fe particles 



SED of the SMC 



Figured] shows the observed global SED of the SMC, after removal of smooth foregrounds 
and backgrounds, from the following compilations: 



Haynes et al.l (ll99ll ) reported global flux densities measured with the Parkes 64m tele- 



scope at 1.4, 2.45 GHz , 4.75 , and 8.55GHz; the 1.4 G Hz flux density is a revision of the 
result of lLoiseau et al.l ( 119871 ). iMountfort et al.l (119871 ) measured the 2.3 GHz flux density 
with the Hartebeesthoeck 26m telescope. 



The TopHat balloon experiment (lAguirre et al.ll2003l ) measured the flux in 4 bands (245 
630 GHz) in a 2.40° radius region (Q = 0.00544 sr) centered on the SMC. Foreg round removal 
was d one by subtracting the mean brightness of adjacent off- source regions. lAguirre et al. 



(120031 ) also extracted 100, 140, and 240/xm fluxes for COBE-DIRBE flSilverberg et al.lll993l ). 



Israel et al.l (120101 ) extracted fluxes for a 2.40° radius region (Q = 0.00544 sr) centered 



on the SMC. We show their extractions for the 10 COBE-DIRBE bands (1.27/im to 248/xm). 
Gordon et al.l (120111 ) extracted fluxes for a 2.25° radius region (Q = 0.00484 sr) centered 
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Fig. 1. — (a) The observed SED of the SMC (2.38° radius region), after subtraction of Galactic foreground 
emission. Estimated contributions from SMC synchrotron and free-free, and from CMB fluctuations are 
shown, (b) SED of the SMC after subtraction of ACMB, free-free, and synchrotron. Flux densities mea- 
sured by IRAS CO B E-DIRBE, TbpHa t, IRAC, MIPS WMA P, a nd Planck are taken from tabulation s by 
Aguirre et all |2003l ). Esrael et al.l poioh . bordon et all l|201ll) . and lPlanck Collaboration et al.l (fcoilbl ) (sec 
text). 



on the SMC, measured using the IRA C (IFazio et al 



cameras on the Spitzer Space Telescope (j Werner et al. 



2004h and MIPS faieke et al.l l2004h 



20041 ). Foreground removal consisted 



of subtracting the mean brightness of an annulus extending from 2.3° to 2.5 C 



Planck Collaboration et al.l ( 1201 lbl ) extracted fluxes for a 2.38° radius region (Q 



0.00542 sr) centered on the SMC. Foreground subtraction consisted of subtracting the mean 
brightness of a 1° annulus aro und the extraction region. We inc lude their extractions for 
Planck 
94 GHz 



(9 bands, 30-858 G Hz IPlanck Collaboration et al.l l2011ah. WMAP (5 bands, 23- 



Bennett et al.l 120031) . and IRAS (4 bands, 12-100/mi; iMiville-Deschenes fc Lagache 



20051 ). IPlanck Collaboration et al.l (j2011b[ ) further corrected the foreground removal by tak- 
ing into consideration the difference in iV(H I) at Galactic radial velocities between the 
background annulus and the extraction aperture^ 

Figure shows the spectrum of the 0.00542 sr region centered on the SMC. We assume 
that the differences in coverage (fi ranging from 0.004 84 sr to 0.00544 sr) are unim p ortant , 
as most of the flux will come from the central regions. IPlanck Collaboration et al.l (|2011bl ) 
estimate that CMB fluctuations add emission corresponding to a mean CMB temperature 



2 We d o not show the "corrected" Planck fluxes for IRAS12 and IRAS25 because the entries for 7^ ub in 
Table 2 of Planck Collaboration et al. ( 2011bl ) do not appear to be correct. 
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excess (ATcmb) = 58yuK over the 0.00542 sr extraction region (relative to the background 
annulus). The spectrum of this CMB excess 



(ACMB)„ = ^(ATcmb) 



dB h 



dT 



(3) 



T=2.726K 



is plotted in Figure [JJa,. 



To isolate the emission from the dust, it is necessary to subtract free-free and synchrotron 
emission. We find the observations to be consistent with synchrotron and free-free spectra 



^isynch 



36. 



11.0 



-1.0 



Jy 



GHz 
ff (lOGHz,T)' 



-h(u-10GHz)/kT 



Jy 



(4) 
(5) 



wi th T = 10 4 K. Our estimate for (10 GHz) is intermediate between the 13 . 4 Jy es timate 



of llsrael et al.l (120101 1 and the 9.05 Jy estimate of iPlanck Collaboration et al.l (l2011bl ). Our 
estimates for Ff? and F^ ndl are shown in Figure [jja. The Gaunt factor gg(u,T) is obtained 
from eq. (10.9) of lDraind (1201 ll ). For n(He + )/n(H. + ) = 1.08 and T = 10 4 K, this corresponds 



to J n e n(B. + )dV = 1.45 x 10 64 cm 3 and an H photoionization rate Qq 



3.7 x 10 51 s _1 . 



The residual after subtraction of (ACMB)„, Fj, and F^ ynch is shown in Figure [Tb. This 
residual is presumed to be emission from dust and (at short wavelengths) stars. A smooth 
curve has been drawn through the points to guide the eye. Subtracting an estimate for the 
starlight continuum as in Figure lb, the integrated A > 5yum dust luminosity of the SMC is 
L d (\ > 5/jm) = 1.00 x lO 8 (£>/62kpc) 2 L . 



Conventional Dust Models 



The observed infrared and submm emission from normal-met allicity star-forming spiral 
galaxies appears to be consistent with physical dust models that were developed to repro- 
duce the observed properties of dust in the diffuse ISM of th e local Milky Way, including 



wavelength-dependent extinction and infrared emission fe.g.. IWeingartner &: Draine 



2001 



Li fc DraineboOll ; IZubko et alfeooi braine fc LfeoO'H : braine fc Fraisse!l2009l : ICompiegne et al 



20101 1 . The models of iDraine fc Lil (120071 . henceforth DL07) consist of amorphous silicate 
grains plus carbonaceous grains; the carbonaceous grains have the physical properties of 
graphite when large, and the properties of polycyclic aromatic hydrocarbons (PAHs) when 
very small. Thi s dust model was a ble to reproduce the global SEDs of the galaxies in the 
SINGS sample (IDraine et al.l 120071 ) . including 17 galaxies with 850/im SCUBA photome- 
try. More recently, the same model has been found to be consistent with both global and 
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spatially-resol ved SEDs of normal- metallicity (i.e., 0.5 < Z/Z & < 2) galaxies in the KING- 
FISH sample f lAniano et aD l2012al jbh . including photometry out to 500/im. 



The DL07 dust mode 



is able to reproduce the observed emission from dust in the 



diffuse IS M of the Galaxy (|Finkbeiner et al.lll999l ) out to wavelengths as long as 2 mm (see 



Fig. 14a of lDraine fc Lill2007l ). However, a significant emission excess (relative to the model) 
appears at A > 3 mm [y < 100 GHz); t his "anomalous microwave emissio n" (AME) has 
been confirmed by numerous observations ( jPlanck Collaboration et al.ll201ld . and references 
therein). T he AME in the Galaxy is th ought to be mainly rotational emission from the PAH 
population ( iDraine Sz Lazarianlll998aU bl). 



DL07 propose that the 3/im < A < 3 mm SEDs of entire galaxies, or large regions within 
a galaxy, can be fit using a dust model consisting of amorphous silicates, graphitic grains, 
and PAHs, and assuming that the dust heating rate is distributed according to 



dMd 
dU 



1 - i)M d , tot 5{U - U miQ ) + 7 M di 



[a -l^- 



iot 



u 



l-a 



ut 



for U m]n <U<U n 



(6) 



where U is the ratio of the local dust heating rate to the heating rate produced by the 
solar-neighborhood starlight radiation field, M d (U) is the mass of amorphous silicate plus 
carbonaceous dust with heating rates < U, and M d)to t is the total dust mass. Eq. (jSJ) is a very 
simple distribution function, with only 5 parameters (M djtot , U m - m , U ma _ x , a, 7), but studies 
of dust emission using this distribution function for the grain heating have been successful 
in reproducing the global e mission from galaxie s, as well as the emission from ~500 pc 
regions within galaxies (e.g.. lAniano et al.ll2012ai ) . The DL07 models are also characterized 
by the PAH abundance parameter ^pah = the fraction of the dust mass contributed by PAH 
particles containing < 10 3 C atoms. 

We vary 5 parameters - the total dust mass M djtot , the PAH abundance parameter (?pah> 
and the starlight heating parameters U m m, a, and 7; t/ max = 10 7 is kept fixed. Because we 
do not include a realistic model for the starlight contribution to the SED, the model is fit 
only to data at A > 3/im, where reddening by dust should be minimal. Because the emission 
at A > 3 mm [v < 100 GHz) may include a substantial contribution from "spinning dust" , 
the DL07 model + starlight is fit only to A < 2 mm data. 

If we allow Um\ n to be as low as 0.2, we obtain Model 1, shown in Figure [2j This model 
has a total dust mass M dtot = 1.3 x 10 6 M Q , exceeding the upper limit of 1.1 x 1O 6 M (see 
Tabled]). Despite using more dust than is allowed, Model 1 provides insufficient emission at 
A > 2 mm. 

Because Model 1 violates the dust abundance limit, we try fitting the DL07 model to 
the same data, but now limiting f/ m ; n > 0.4. The resulting Model 2 has a total dust mass 
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\(/j,m) A((im) 

Fig. 2. — Data: Observed SED of the SMC, after removal of CMB fluctuations, and subtraction of free-free 
and synchrotron emission (see Figure [TJ. Solid line: model consisting of the sum of T = 5000 K starlight 
(dotted line), emission from the DL07 dust model (dot-dash line), and a spinning dust component peaking 
at 40 GHz (dashed line), (a) The total dust mass exceeds the maximum allowed by a factor ^20% (sec Table 
[I}, (b) A model with the dust mass within allowed limits. In both models, the spinning dust component has 
been adjusted to reproduce the observed 20-50 GHz emission. Both models provide insufficient 60-300 GHz 
emission. The models in (a) and (b) have dust luminosities L d = 7.0 x 1O 8 L and 6.9 x 10 8 L Q . 



that does not violate the upper limit in Table (TJ but the quality of the fit to the SED is 
somewhat worse than for Model 1, with an even larger deficiency at A > 2 mm (see Fig. [2b). 



5. Spinning Dust 



We can add a spinning dust component to raise the emission in the 20-60 GHz range. 
The spinning dust emission in the diffuse ISM of the Galaxy peaks at about 40 GHz; we 
expect a similar peak frequency in the SMC. What level of emission is expected for spinning 
dust in the SMC? 



Draine fc Lazarian ( 1998al lbl) argued that the anomalous microwave emission in the 



,•(*!) 



(40GHz)/n H 



MW 



lxlO- 17 Jysr^cnrTf- 1 , 



Galaxy, with an observed emissivity per H nucleon 

is primarily rotational emisssion from the PAH population. The PAH abund ance is mea- 
sured by (/pah- Dust in the solar neighborhood is thought to have (/pah ~ 4.7%. iLi fc Draine 
( 120021 ) found that (/ pah in the SMC was spa tially variable and, on average, much lower than 
in the Milky Way. ISandstrom et al.l (120 lOf ) confirmed this, finding a mean ((/pah) ~ 0.6%. 
We expect the dust/gas ratio in the SMC to be lower by about a factor ~Zsmc/-Z© ~ 0.25. 
Therefore the PAH abundance per H is down by about a factor ~(0.6/4.6) x 0.25 = 0.033. 



9 



Thus we estimate the spinning dust emission in the SMC to be 
~jt d) (40GHz " 



n H 

J SMC 

AFj sd) (40GHz) 



0.033 x 1 x 10" 17 Jysr^cn^H" 1 « 3.3 x 10~ 19 Jy sr" 1 cm 2 H~f7) 
jt d) (40GHz) 



n H 



J SMC 



M H n 

x— ZT 2 ^5Jy . 
m H 



(8) 



In Figure [2] we have added an emission component with a spectrum! 3 ] 

AF (sd) = AF (sd )(z/o) exp ^ _ (1//vb) 2] 



(9) 



representative of what is expected for spinning dust. If we set uq = 40 GHz and AFj sd) (40 GHz) 
5 Jy - consistent with the estimate in Eq. (jSJ) - the 20-50 GHz observations are accounted 
for, as seen in Figure [2j 

While spinning dust appears able to account for the observed 20-50 GHz emissi o n, the 
emission between 50 and 300 GHz remains much stronger than expected. iBot et al.l (120101 ) 
suggest that the 50-300 GHz excess could also be due to spin ning dust emission. How- 
ever, theoretical models of rotational emiss i on from small grains (|Draine fc Lazarianl Il998bl ; 
Ali-Haimoud et alJliooa iHoang et aflboiol boill ; Isilsbee et aliboilh predict little rotational 
emission above ~100GHz unless hi gh densities and wa r m gas temperatures are present in 
the emitting regions. For example, iDraine fc Lazarianl (jl998bl . see Fig. 13) calculated the 
spinning dust emission from a model PDR with tt-h = 10 5 cm -3 , T = 300 K, illuminated by a 
radiation field U ~ 3000. The PDR was assumed to have abundances of small grains relative 
to big grains reduced by a factor of 5 relative to diffuse clouds in the solar neighborhood, 
approximating the observed reduction in (/pah in the SMC. Viewed face-on, the total IR lumi- 
nosity/area m 3.6 x 10~ 3 £7 ergscm~ 2 s _1 = 11 ergs cm" 2 s^ 1 . The spinning dust emis- 
sion for this model peaked near 110 GHz, with (z/L1 s< ^)iooghzA4 



Thus Ltih/(uL 



(sd) . 

100 GHz ~ 



7.4 x 10 7 ergs cm 2 s 1 . 



1.5 x 10 
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At 100 GHz, the model in Figure Wp has a deficit AF„ 25 Jy, corresponding to 
iooghz = 47rD 2 (^AFy) 100 GHz = 3OOOL . To account for this would require PDRs 
with a luminosity Ljr = 4.5 x 10 10 L o - completely inconsistent with the observed Ltir = 
1 x 10 8 L o . It is evident that spinning dust cannot account for the observed 50-200 GHz 
emission in the SMC. Here we consider magnetic dust grains as an alternative. 



This simple f orm, adequate for the present purposes, approximates the spectra obtained by detailed 



calculations (e.g., Draine fc Lazarian 1998b : Ali-Hai'moud et al. 2009t Hoang et al. 2010l 2011 : Silsbee et al 



2011 



Ysard et al.ll201lh . 
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A( J um) X(ixm) 

Fig. 3. — Similar to Figure [21 but with metallic Fe nanoparticles added to the dust model. The Fe particles 
are assumed to be at T = 40 K in Model 3 (panel a) and T = 20 K in Model 4 (panel b) . 
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Fig. 4. — Similar to Fig. [3] but with emission from (a) maghcmite (Model 5) and (b) magnetite (Model 6) 
nanoparticles at T = 17 K. 

6. SMC Dust Models Including Magnetic Dust 

Because magnetic materials have enhanced absorption at microwave and submm fre- 
quencies, it is of interest to see whether the mm- and cm-excess seen in the SMC could be 
due in part to thermal emission from magnetic grain materials. In Figures [3] and H] we model 
the observed emission from the SMC as the sum of three components: "normal" dust (the 
amorphous silicate, graphite, and PAH model of DL07), a population of magnetic grains, and 
spinning dust. In each case, the spinning dust contribution is assumed to peak at 40 GHz, 
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with the peak flux density adjusted to fit the observations in Figure [T^,, giving a reasonably 
good fit in the 20 — 60 GHz region. 

Metallic iron nanoparticles are introduced in Fig. [3J We consider Fe grain temperatures 
of 40 K (Fig. (3^l) and 20 K (Fig. [3b). If the Fe nanoparticles are, fo r the most part, free-flier s 



heated by typical starlight, then T « 40 K is expected (see Fig. 4 of lDraine fe Hensleyll2012l ). 
If, on the other hand, the Fe nanoparticles are inclusions in larger composite grains, then the 
T ~ 20 K temperature is appropriate, consistent with the temperature of the "normal" dust. 
In each case, the Fe grain abundance is adjusted to reproduce most of the observed emission 
near 100 GHz, then a model using DL07 dust is used to provide the additional emission 
required to reproduce the observed SED at shorter wavelengths, and finally a spinning dust 
component peaking at 40 GHz is added to bring the model into agreement with the 20- 
50 GHz observations. 

We also consider nanoparticles of maghemite (Fig. Hk) and magnetite (Fig. |4b). For 
these we assumed temperatures T ps 17 K consistent with being inclusions within nonmag- 
netic dust grains. 

The model with mag hemite (Fig. W) has M Fe = 2.2 x 1O 5 M of Fe in maghemite (total 
maghemite mass 3.1 x 1O 5 M ) and the model with magnetite (Fig.Hb) has Mp e = 2.2xlO 5 M 
of Fe in magnetite (total magnetite mass 3.0 x 1O 5 M ). The Fe mass, and total dust mass, 
does not violate the mass budget (see Tabled]). We conclude that the observed mm- wave 
emission from the SMC can be accounted for by models with reasonable abundances of 
normal dust plus metallic Fe, maghemite, magnetite, or some combination of these three 
materials. 

If the nanoparticles are present as inclusions in larger grains, it is clear that the size 
distribution of the larger particles can be adjusted to be compatible with the observed 
wavelength-dependent extinction in the SMC. But is it possible for the bulk of the interstellar 
Fe to be in free-flying nanoparticles? We have calculated the extinction contribution in the 
optical and UV, assuming that 100% of the Fe is in particles of a singl e type , and using 
dielectric functions for Fe, Fea04, and 7-Fe203 from iDraine fe Hensleyl (J2012J). Figure 



show s the calculated ext inction per H, together with the observed extinction in the SMC 



Bar ([Gordon et al.ll2003l ). In no case does the calculated extinction/H exceed the observed 
extinction. Therefore, the observed extinction is not incompatible with the possibility that 
much of the Fe is in free-flying nanoparticles. 

Nevertheless, we consider it most likely that the bulk of the magnetic nanoparticles 
would be present as inclusions in larger grains, since we know that most of the grain mass 
is in grains with radii a > 0.1/xm 
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Fig. 5. — Extinction contributed by nanoparticles if 100% of the Fe is in particles of Fe, Fe3C>4, or 7- 
Fe2Q3, and radii a = 5 nm or a = 20 nm. Also shown (symbols) is the observed extinction in the SMC Bar 



(jGordon et al.l 120031) . The observed extinction does not rule out the hypothesis that most of the Fe is in 



free-flying nanoparticles. 



Discussion 



7.1. Solid-Phase Iron in Low Metallicity Galaxies 



At submm - mm frequencies, the SED of the SMC is significantly less steep than the 
SED of normal-metallicity spiral galaxies, including the Galaxy. If the enhanced emission 
of the SMC at mm-wavelengths is due to Fe or Fe oxide grains, then such grains must 
account for a larger fraction of the dust mass in the SMC than in normal-metallicity spirals: 
conditions in the SMC must be in some way more favorable for their production or survival 
than in normal star-forming galaxies. 

Spitzer Space Telescope observations of globular clusters have detected excess infrared 
emission from the most luminous giant stars, indicative of dusty winds. In many cases, the IR 
spectrum of the infrared excess is dominated by a featureless continuum at A > 8/xm. Glob 



ular cl usters wher e such featureless s pectra have been s een include 47 Tuc f|McDonald et al. 
2011ah . NGC 362 faover et al.ll2009h . and Omega Cen jMcDonald et al.l boilbh: these three 



clusters have metallicities [Fe/H]» —0.7, —1.3, and —1.5, respectively (jHarrisI Il996l . 2010 



edition)0. The featureless emission might be attributed to hot amorphous carbon grains, 
but carbon solids are not expected to form in these oxygen-rich outflows. The featureless 
continuum has therefore instead been attributed to thermal emission from metallic Fe grains 



4 The updated Harris cluster catalogue can be found at 



http : //www . physics . monaster . ca/ $\sim$harr is . mwgc . dat 
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([McDonald et al.ll2010U2011al Jbl). Thus, low-metallicity AGB stars provide a possible source 
for metallic iron or Fe oxide grains in low-metallicity galaxies such as the SMC. 

T he iron-rich ejecta of Type la supernovae constitute a second potential source of iron 
grains (jDwek!ll998l ). However, to date there is no evidence of d ust formation in SN la ejecta, 
despite sensitive searches toward the Tycho and Kepler SNRs (IGomez et al.ll2012l ). 



Type II supernovae are known to form dust in the ejecta in at least some cases (ISugerman et al 



20061 ; iMatsuura et al.ll201ll ). and i t is conceivable t hat Fe-rich portions of the ejecta might 
condense metallic Fe or Fe oxides. iRho et al.l ( 120081 ) made models to reproduce the 5-38/im 
spectra of the Cas A ejecta; their global model had O.O28M of dust, of which 37% was 
metallic Fe. 



Baron et al.l (119771 ) observed that lunar soil grains have an increase in the concentration 
of Fe near the surface, with some of the Fe in metallic form. These surface layers ("rims") 
reflect exposure of the grains to cosmic rays and the solar wind. "Inclusion-rich rims" consist 
of an amo rphous silica-rich matri x with abundant metallic Fe inclusions, typically < 10 nm in 
diameter (jKeller fe McKay! 1 19971 ). Inclusion- rich rims are compositionally distinct from the 
host grain, and are thought to have formed by deposition of atoms from vapors produced by 
nearby sputtering or impact events. In the laboratory, irradiation of olivine by 4 keV He ions 
is obse rved to lead to alt e ration of the surface layers, with forma tion of metallic Fe nanopar- 
ticles (IDukes et al.lll999l ; ICarrez et al.ll2002t iLoeffler et al.ll2009l ). Metallic Fe nanoparticles 
are found as inclusions in inte rplanetary du st particles known as GEMS ( "Glasses with Em- 
bedded Metals and Sulfides"; lBradleyill994l ) . Thus, it is reasonable to consider that some of 
the Fe in interstellar grains may be in metallic Fe inclusions. 

Fe-rich grain material is injected into the ISM from stellar sources; as seen above, 
this may include metallic Fe. Additional conversion of gas-phase Fe to solids must take 
place in the ISM to account for observed low gas-phase abundance of Fe, particularly in 



(Draine & Salpeter 


1979a: 


Jones et al. 


1994; 


Draine 


2009) 



dust" (material condensed in stellar outflows) accou nts for only a small fraction - perhaps 



10% - of the interstellar grain mass in the Galaxy (jDraind Il990l . 120091 ). and this is likely 
the case for all galaxies where a substantial fraction of the refractory elements (Mg, Si, Fe) 
is in grains. In such galaxies, including the SMC, the bulk of the grain material must have 
undergone conversion from gas to solid in the ISM. The character of the interstellar dust 
will therefore be largely determined by interstellar processing. 

Sputtering by energetic H and He can alter the composition of interstellar dust . Sput- 



tering yields have been discussed by a number of authors (e.g., iDraine fe Salpeterl Il979bl ; 
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Tielens et al.lll994f ). For a composite material, sputtering yields for H and He will be larger 
for the lighter elements in the target, and sputtering will therefore leave the surface layers 
enriched in heavy elements (such as Fe). The grain material that survives sputtering will 
therefore become Fe-rich, perhaps even metallic Fe. Studies of element al depletions in the so- 



1997; 


Jenkins 


2009) 


Weltv et al. 


(2001) 



Jenki ns 2009). Based on the observed depletion patter toward Sk 155 in the SMC, 



sightline to Sk 155) is in the form of metallic Fe or Fe oxides. 

Rates for grain growth by accretion are proportional to the metallicity, while rates 
for grain destruction by H and He sputtering are not. The balance between grain growth 
and destruction, and the composition of the extant material, will therefore depend on the 
metallicity of the ISM. Hence, the apparent difference in grain composition between normal- 
metallicity spirals (like the Milky Way) and low-metallicity dwarf galaxies such as the SMC. 



7.2. High-Frequency Magnetism and the Gilbert Equation 

The models pr esented here use absorpt ion cross sections C a b s (w) for m agnetic grains 



calculated following iDraine fe Hensleyi (120121 ) . who used the Gilbert equation ( lGilbertll2004l ) 
to model the frequency-dependent magnetic response of Fe, maghemite, and magnetite. The 
Gilbert equation uses an adjustable dimensionless parameter oq to characterize the dissipa- 
tion. We have adopted «g ~ 0.2 for purposes of discussion, but the existing experimental 
literature employs a range of values of oq. If ac were to be smaller than 0.2, the v > 100 GHz 
opacity would be reduced, and the mass of Fe required to reproduce the observed emission 
of the SMC would correspondingly increase. If olq J$ 0.05, magnetic grain models to explain 
the ~3 mm emission would be ruled out by abundance constraints. 

Quite aside from the question of what value to use for qq, it is also important to rec- 
ognize that the prescription for dissipation in the Gilbert equation, while mathematically 
convenient, is not based on an underlying physical model. Empirical evidence for the ac- 
curacy of the Gilbert equation at high frequencies is scant. Laboratory measurements of 
electromagnetic absorption in Fe and Fe oxide nanoparticles at frequencies up to 500 GHz 
are needed to validate use of the Gilbert equation at these frequencies. 



15 



7.3. Polarization 

Based on starlight polarization studies, the magnetic field in the SMC appears to lie pr i- 
marily in the plane of the sky, with substantial large-scale coherence (IMao et al.ll2008l . 120121 ) . 
While our understanding of the physics of grain alignment remains incomplete, dust grains 
in the SMC are expected to be partially aligned with long axes tending to be perpendicular 
to the local magnetic field Bo- Electric-dipole emission, which dominates for A < 500/im, 
will be polarized with E w _L B . 

As seen above, magnetic dipole emission from magnetic nanoparticles may become im- 
po rtant for A > 1 mm . The polarization of the magnetic dipole emission has been discussed 
by iDraine fe Hensleyl (|2012| ). If the magnetic nanoparticles are free- fliers, and are aligned 
by a Davis-Greenstein-like mechanism, the magnetic dipole emission will be polarized in the 
same sense as the FIR emission, b ut the fractional po l arizat ion may be even larger than that 
of the FIR emission (see Fig. 9 of IDraine fe Hensleyl ( 120121 )). Alternatively, if the magnetic 
nanoparticles are present as randomly-oriented inclusions within larger aligned grains, the 
magnetic dipole emission will be polarized with E w \\ Bq. As a result the fractional polariza- 
tion may decrease by a factor ~2 as the frequency decreases from 200 GHz to ~40 GHz, and 
for magnetite, maghemite, or Fe spheroids, the net polarization undergoes a reversal (i .e, 
changes from E w _L B to E w || B ) near ~15GHz (see Fig. 10 of IDraine fe Hensleyl 120121 ). 



PlanckwiW measure the polarization at 30, 44, 70, 143, 217, and 353 GHz. Unfortunately, 
there are two additional factors that will complicate interpretation of the dependence of 
polarization fraction on frequency: 



Emission from spinning dust becomes increasingly important with decreasing fre- 
quency, peak ing near ^40 GHz. If th is emission component is minimally polarized, 



as predicted (ILazarian fe Drainell2000l ). it will cause the fractional polarization to de- 
crease with decreasing frequency in the 60-100 GHz region. 

There may be more than one grain type contributing to the normal "electric dipole" 
emission at lo ng wavelengths, as in t he mixtures of silicate and carbonaceous grains 
considered by IDraine fe Fraissd (120091 ). In this case, even the "normal" electric dipole 
emission alone may have the fractional po l arizat ion depending significantly on fre- 
quency. In the models of IDraine fe Fraissd (120091 ). the fractional polarization is pre- 
dicted to increase with decreasing frequency. 



Actual reversal of the polarization below ~15 GHz would be an unambigous indication of 
magnetic dipole emission from magnetic inclusions, but this may be overwhelmed by the 
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increasing importance of synchrotron emission (polarized with E w _L Bo) as the frequency 
falls below 10 GHz. 



8. Summary 

The principal conclusions of this paper are as follows: 

. We show (see Figs. E] and SJ) that the SED of the SMC can be approximately reproduced 
by a mixture of "normal" dust (illuminated by a plausible range of radiation intensities) 
plus emission from a population of small (a < O.Olyum) magnetic nanoparticles. We 
consider three magnetic materials: metallic Fe, magnetite Fea04, and maghemite 7- 
Fe203. It appears that any of these 3 materials, or a combination of them, can provide 
enough emission at A > 1 mm so that a combination of "normal dust", spinning dust, 
and magnetic dust can acount for the observed SED of the SMC. 

. If conditions in the SMC are conducive to a large fraction of the interstellar Fe being in 
magnetic nanoparticles, other low metallicity galaxies may also have mm-wave emission 
dominated by magnetic dipole emission. 

. While it seems natural for the magnetic nanoparticles to be inclusions in larger grains, 
the observed extinction does not rule out the possibility that the magnetic nanoparticles 
might be independent free-fliers. 

. If the magnetic nanoparticles are present as randomly-oriented inclusions in larger 
silicate grains, the polarization is expected to fall as the frequency decreases below 
~ 200 GHz. It may be possible to test this prediction with measurements by Planck of 
the polarized emission from the SMC. 



Our m odels are based on high-frequency magnetic properties as estimated by lDraine fe Hensley 



( 120121 ) using the Gilbert equation. Laboratory studies of the high-frequency {y > 
100 GHz) magnetic properties of metallic Fe and Fe oxides are needed to improve our 
understanding of magnetism at high frequencies. 
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